Separation of variables on a non-hyperelliptic curve 



V.G. Marikhin and V.V. Sokolov 
: February 8, 2008 

O 

o 

D 

Q ■ 1 Introduction 

r-> 

In the paper we consider several models that admit a separation of variables on the following 
. algebraic curve of genus 4: 

^ ' m Y) = s 6 Y* + /(OF 4 + k(0Y 7 - S(0 = 0, (1.1) 

where 

! ^6 = S ^S VI (0, HO = ^"(0+4K 2 +e 2 ?+e 1 ), 1(0 = ^S IV (0-jk"(0- (1-2) 

^ ■ Here a, 5 are parameters, S is arbitrary sixth degree polynomial. If 5 = 0, then the curve is 
^ \ hyperelliptic of genus 3. 

The class of curves (jl.lj) with 5^0 can be described as follows. Let P(?7,£) = be an 
£^ ■ arbitrary cubic. In the generic case it is an elliptic curve. Then (jl.lj) is a double cover over 
^ the cubic defined by the formula r\ = £ 2 — 5 2 Y 2 . 

X : 

5_| ■ To present out approach, we start with the Hamiltonian 

H = apf + cp\ + dpi + ep 2 + /, (1.3) 

where 

4s 2 S{ Sl ) 4si5(s 2 ) j J J 
a = , c = , a = —si , e = — s 2 - 



S\ — S 2 Si — S2 Si — S2 Si — ^2 

£ 2 s 2 S"( Sl ) - Sl S"(s 2 ) _ 6 2 s 2 S'(s 1 ) + s 1 S'(s 2 ) | 
40 si - s 2 4 (si - s 2 ) 2 

35 2 s 2 S( Sl ) - Sl S(s 2 ) 



;i.4) 



(si - * 2 ) 3 



+ asis 2 , 



1 



and 

si - s 2 

as the most general model. It can be easily verified that H commutes with the function 

„2 , ri2 



where 



K = Api + Cp 2 2 + D Pl + Ep 2 + F, (1.5) 



A= 4sw C = _^M D= ' E= _L 



S\ — S 2 S\ — S2 S\ — S2 S\ — S2 

P S»{s 2 )-S»{s 1 ) 6* S'( Sl ) + S'(3 2 ) 35* S(s 2 ) - S( Sl ) 

F = J7{ + T — 7 ^— + ~a 7 V3 a>(s 1 + s 2 ) (1.6) 

40 Si - s 2 4 (si - S2) 1 4 (si - 52)^ 

with respect to the standard Poisson bracket {p a , sp} = 5 a p. If 5 = 0, then the terms in H 

and K linear in momenta disappear and H belongs to the Stackel class. In this case si,s 2 

are separated variables. The corresponding Abel's transformation on the hyperelliptic curve 

of genus 3 has the form 

dsi ds 2 , Sidsi s 2 ds 2 

+ ; = dt, , + ; = 0. 



Here 

T(s) = 16S(s) (as 2 + e 2 s + ei), 

and the constants ei and e2 are values of the integrals H and K, respectively. 

The main result of the paper is an explicit formula for the action in the case 5^0 
depending on two parameters. We find it directly from the Hamilton- Jacobi equation. 

In Sections 4-6 we find the action and a separation of variables for the Kowalewski 
hyrostat, the Clebsch and the so(4) Schottky-Manakov spinning tops in a way analogous 
to that used in Sections 2-3 for the Hamiltonian ()1.3)) . Many significant papers (see, for 
instance, [1111101211311111131111111112111301211 and references there-in) are devoted to 
these models. In particular, it is known (2131211) that the models are linked via some changes 
of variables. However according to a common opinion of experts, a "satisfactory" separation 
of variables for the above models is still not found. 

For the Clebsch and the so(4) Schottky-Manakov tops S is a polynomial of degree 4. 
An appearance of the corresponding algebraic curve of genus 3 could be predicted: the 
characteristic curve for the 3 x 3-Lax operator found by A. Perelomov (3] for the Clebsch 
top has the same form. However in this paper we don't use any Lax representations and the 
curve arises during computations in a rather natural manner. 

The case S = const also is of interest. The inhomogeneous hydrodynamic type system 
(12. 7|) corresponding to the pair (jl.3|) . (jl.5j) coincides with the Gibbons-Tsarev equation 0. 
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The separation of variables found in this paper gives rise to a family of elliptic solutions for 
the Gibbons- Tsarev equation (see Section 7). 

Acknowledgments. The authors are grateful to A.V. Borisov, Yu. N. Fedorov, E.V. 
Ferapontov, V.B. Kuznetsov, and F. Magri for useful discussions. The second author (V.S.) 
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2 Equations of motion 

The Hamilton equations corresponding to (jl.3|) and (jl.5|) have the form 

— - = 2a Pl + d, -f = 2cp 2 + e. 2.1 

at at 

and 

^ = 2A Pl + D, ^ = 2Cp 2 + E. (2.2) 

dr dr 

Finding p±,p 2 from (|2.1j) and substituting them in the relations H = e\ and K = e 2 , we get 

<?' 2 s 2 4^ 2 
(*i " ^[si^ - ^2^] + 16( ei - f)s lS2 + -(si S( Sl ) - si S(s 2 )) = (2.3) 

0[Si) J{ s 2) {Si — S 2 ) 

and 

AX 2 

lGs 1 s 2 ((e 1 -f)(s 1 +s 2 ) + (e 2 -F)s 1 s 2 )—. ^ (si S( Sl )( Sl - 2s 2 ) - s\ S{s 2 ){s 2 - 2s x )) = 0, 

C»i - s 2 y 

where / and F are defined by ()1.4j) . f)1.6|) . 

One of the basic technical problems is a rewriting of this system in a compact form. The 
straightforward elimination say s 2 leads to a cumbersome equation of degree 4 with respect 
to s\. 

Proposition 1. Let (u, v) be any solution of the following system of equations: 

L(s 2 )u 2 -2L{ Sl )v- M{s 2 , Sl ) = 0, L{ Sl )v 2 -2L{s 2 )u-M{ Sl ,s 2 ) = 0, (2.5) 

where 

M(x, y) = 3L(y) + L'(y)(x - y) + k(y)(x - y)\ L(x) = 5 2 S(x). 
Then the derivatives 

s 2 u + s 1 . T s 1 v + s 2 

si = -J , s 2 = -J (2.6) 

Si - S 2 Si - s 2 
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satisfy the system ()2.3|) . (|2.4|) . 

The elimination of variables pi,p2 from (J2.1)) . (|2.2jl leads to the inhomogeneous hydro- 
dynamic type system 



{si)t + s 2 {s 1 ) T = -J, (s 2 ) t + s 1 (s 2 ) T = J. (2.7) 



It follows from this formula and from (J2.6j) that 



( Sl )r = J^~, (S 2 ) T = J^l. (2.8) 
Si -s 2 Si - s 2 

Moreover, we find from (|2.6|) that 

Ju Jv 
ob(si) ob(s 2 ) 

The following parameterization for solutions of ()2.5|) turns out to be crucial for the 
separation of variables. 

Proposition 2. Let z\, z 2 , z% be solutions of the equation 

z 6 L(s 2 ) 2 - L(s 2 )M(s 2 , Sl )z 4 + M( Sl ,s 2 )L( Sl )z 2 - L{ Sl f = (2.10) 

such that 



L(x) M(y,x) _ 2 , 2 M{x,y) 111 



Then 
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1, 


— + 


— + 
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1 


^2 
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1, 
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Z3 


Z\ 


2^2 



(2.11) 



are so 



lutions of (12.5)1 . 



3 Separation of variables. 

To separate the variables in the case 5 ^ 0, we find the action function S{s\, s 2 ) in an explicit 
form. 

Let us consider the system of equations 

*(£,¥) = 0, F 2 =i(£- Sl )(£-s 2 ), 



where the polynomial Y) is given by (jl.lj) . It is easy to verify that if we substitute to 
the first equation the expression for Y 2 taken from the second equation, then leading powers 
of £ cancel and the equation for £ turns out to be cubic. By £i(si,s 2 ), i = 1,2,3 denote 
the roots of this equation. 

Theorem. The function 



1 3 

S(s 1 ,s 2 ,e 1 ,e 2 ) = 



4 

n=l 



5n 

X + u £n ~ 2 (Sl + S 2 ) f d£ 

o arctanh ■ 



5Y(£ n ) J y(0 



(3-1) 



satisfies the Hamilton- Jacobi equation 



dS dS \ dS dS 

where the functions H(pi,p2, si, S2) and K(pi,p 2 , si, s 2 ) are defined by formulas (jl.3J) - (J1.5J) . 
Proof. It follows from (|2.9|) that the partial derivatives of the action S are given by: 

dS J u dS J v 

Ihl = 85(si)' dT 2 = ~8S{s 2 )' 

The compatibility condition for this system has the form 

S(s 2 )—{Ju)+S( Sl )—{Jv) = 0. 

OS 2 OS\ 

The latter relation follows from the identities 

{L'{y) + M x {x, y)) {x - y) + 6L{y) - 2M{x, y) = 0, 

(3.3) 

{x - y) (M y {x, y) + M x (y, x)) + 4(M(x, y) - M{y, x)) = 0, 

which can be verified straightforwardly. 

Consider an auxiliary Hamiltonian of the form (compare with formula (|2.1(J)0 

H = p\L(s x ) + p 2 lP2 M(s 2 , Sl ) + Pip 2 2 M( Sl , s 2 ) + p\L(s 2 ). (3.4) 

It can easily be checked that identities ()3.3)1 are fulfilled if the condition {H, K} = 
holds, where 

K = -(si - s 2 ) 2 pip 2 . 
It is easy to verify that if {pi, Sj} = 5ij, then the Poisson brackets between functions 

a =pi+ p 2 , ai = sipx + s 2 p 2 , a 2 = s\p x + s\p 2 (3.5) 



are given by 

{a ,ai} = a , {a ,a 2 } = 2a 1 , {a ll a 2 } = a 2 . (3.6) 

The expression Q = a\ — a^a 2 is a Casimir function for the linear s/(2)-brackets ()3.6j) . It 
turns out that K rewritten in variables (|3.5|) coincides with Q. The Hamiltonian H can also 
be expressed in terms of variables (|3.5|) only: 



H = a L — - S (aa 2 + e 2 a\ + eioo) + , 

\a J 4 a 

It is seen from this formula that the variable 

ao Pi + P2 

should play a key role in a description of properties of polynomials (|2.1()|) . (|3.4jl . Suppose 
H(pi,p 2 ) = 0; then a solution z for equation (j2.10|) can be expressed in terms of £ by the 
formula 

Equation (|2.10|) is cubic with respect to £. It can be rewritten as (jl.lj) . where 

^ 2 = f (3-8) 

Our task is to find the solution of ()3.2|) explicitly. According to formulas (j2.11|) . it suffices 
for that to solve the system of equation 

da J 




dsi 8S( Sl ) 4(si - s 2 

(3.9) 

da J 1 5 

fc~2 ~ ~8S(s 2 )^ ~ ~4(si-s 2 ) 

Here z(si, S2) is arbitrary root of equation (j2.1()j) and £ is the corresponding (see (J3.7j) ) value 
of £(si, s 2 ). The action function 5 required can be obtained as the sum of the three solutions 
for ()3.9|) . corresponding to the three branches of the function £(si, s 2 ). 

First, we find a function ao such that ()3.9|) is fulfilled under condition that £ is a param- 
eter, which does not depend on Si,s 2 . It is not hard to see that 

cr (4) s i) s 2j = 70 arctanh —== 



4 ^_ Sl )(^_ a2 ) 



Taking into account the formula 



d 1 

^7^0 = T 



V(e-^i)(e-^ 



and replacing a/(£ — ~ s 2 ) by we obtain expression ()3.1j) for the action function. 
Differentiating the action with respect to the parameters and taking into account that 

dS _ dS 

o — — t + Ci, — — — C 2 , 

oe\ ae 2 

we finally get 

3 3 

dt = ^w 1 (e„), = J> 2 (£ n ). (3-10) 

n=l n=l 

Here iO\, u; 2 belong to a basis 

wi(0 = — , w 2 (0 = — , w 3 (0 = , w 4 (0 = (3.11) 



def 



of holomorphic differentials on the curve Y) = 0. In ()3.11|) we use the notation Z = |p. 
The differential 004 plays a special role. In variables (£,??), where 

r 1 = e-8 2 Y\ (3.12) 

the curve (jl.ljl becomes the following cubic 

W+s 5 ?? 2 £+^4?7(?7+4£ 2 )+is3^+2^^ 

(3.13) 

The holomorphic differential of the elliptic curve ()3.13j) coincides with U4 up to transforma- 
tion fling). 

It follows from ()3.8j) that the functions are linked by the relation 

1= + + ?m (3 i 4 ) 

S 2 (6 - &)(& - &) (6 - 6) (6 - 6) (6 - &)(& - 6) ' 1 ' ' 

which allows us to determine the two functions Si(£), s 2 (£) starting from £i(t), £2^), ^(t)- 
Thus we have three conditions ()3.1()|) . ()3.14j) to determine three functions £«(£). 

Condition 1)3.14)1 can be rewritten in the variables (£, rf) as 

mfa - 60 + - 6) + - 6) = o. 

This means that the corresponding points r^) belongs to the intersection of elliptic curve 
()3.13jl and the straight line rj = £(si + s 2 ) — sis 2 . 
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4 Schottky-Manakov spinning top. 



It is well-known (see, for example ) that the Hamiltonian 

H = (St, AS,) + 2(S U BS 2 ) + (S 2 , AS 2 ), 

where A = diag(ai, a 2 , a 3 ), B = diag(bi, b 2 ,b 3 ), commutes with a quadratic polynomial K 
with respect to the spin Poisson brackets 

{S?,Sf} = K e af > Y SZ6 ij (4.1) 

iff 

b\(a 2 - a 3 ) + bl(a 3 - at) + b\(ai - a 2 ) + (at - a 2 )(a 2 - a 3 )(a 3 - at) = 0. (4.2) 
Here e a ^ is the totally skew-symmetric tensor, k is a parameter. 

Since H and K may be replaced by arbitrary linear combinations of H, K and the 
Casimir functions 

Ji — (Si, Si), J 2 = (S 2 , S 2 ) 
for the brackets (|4.1j) . the integral K can be reduced ^Hl to the form 

K = 2(Si,CS 2 ), C = diag(ai, a 2 , a 3 ). 

Without loss of generality the matrices A and B defined the Hamiltonian H can be choosen 
as follows 

A = diag(— a\, —a\, — a^), B = diag(a 2 c*3 + A«i, a 3 ai + Xa 2 , aia 2 + \a 3 ), 

The arbitrary parameter A corresponds to the shift of H by A K. 

Reduction to the standard brackets. Let us fix the values of the Casimir functions: 
(Sk, Sk) = jl- Then the formulas 

S k = p k K(q k ) + ^K'(q k ), where K(q) = ((q 2 - 1), i(q 2 + 1), 2q), (4.3) 

define a transformation of the Poisson manifold with coordinates St, S 2 and brackets (j4.1|) . 
where k = —2i, to the manifold with coordinates Pi,p 2 , qt, q 2 and canonical Poisson brackets 
{Pa> qp} = 5 a p. As a result of this transformation we get 

H = pjr(qi) + —pi r'(qi) + —r"(q x ) + p\r(q 2 ) + — p 2 r'(q 2 ) + -|r"(g 2 ) + 

(4 4) 

2 Pi + -W^T~ P2 + Z (Qu 92), 



2 dqij V 2 % 
8 



K = 2 ( Pl + £ L + ^JL ) g2 ), 



2 3gi / \ 2 <9g 2 
where 

r(x) = f (K(x), AK(x)) = -a\(x 2 - l) 2 + a 2 (x 2 + l) 2 - 4a 3 x 2 , 



(a 2 a 3 + Aai)(x 2 - l)(y 2 - 1) - {a 3 ai + \a 2 )(x 2 + l)(y 2 + 1) + 4(aiQ;2 + Aa 3 )xy, 

y) d = (K(x), CK(y)) = a^x 2 - l)(y 2 - 1) - a 2 (x 2 + l)(y 2 + 1) + 4a 3 xy. 
It is possible to check that 

Z 2 (x, y) - r{x)r{y) = W(x, y)W(x, y), 

where W is a polynomial quadratic in each of variables. 
Remark. Note that the more general Hamiltonian 

describing an interaction of iV spins, can be reduced to 

H = J2 9 ij PiPj + Yl aipi + v ' 

where 

a « = ^ir"(x')ir"(^)c«, 



dg 1 . 5(7* 



dx k 2 * dx l ' 

(9 2 ^ ifc 1 . 2 <9 2 # fcfc 



by a similar transformation 



4 ^— ' dx' l dx k 6 ^— ' dx k dx k 

ik k 



S k =p k K(x k ) + 3 ^K'(x k ), k = l,...,N. 
In terms of the coordinates and the velocities the Lagrangian and the energy read as 
L = \ Y ~ ~ aj )- v i E = \ Y 9i0x j - aV) + v, 

ij ij 

where Y.9ij9 jk = 

3 
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An additional integral of the form 

K=j2 Gij PiPi+J2 Ai p* +v 

ij i 

exists iff 

y d& gM _ dg^ Qkl y _ 2 d^ Qkj _ Ak d_f_ + 2 d^ gkj = Q) 

' dx k dx k ' dx k dx k dx k dx k 

k k 

y a kd A. _ A k — + 2g lk — - 2G lk — = 0, V a k — - A k — = 0. 
' <9x fc ftr fc <9x fc ^— ' dx k dx k 

k k 

Diagonalization of the quadratic part. Under transformation ()4.3|) the functions H 
and K take the form 

H — ap\ + 2bpip 2 + cpl + dpi + ep 2 + /, (4.9) 

K = Ap\ + 25pip 2 + Cp* + D Pl + £p 2 + F, (4. 10) 

where the coefficients are some (in our case rational) functions of the variables qi,q 2 - 

General formulas related to a pair of commuting Hamiltonians quadratic in momenta 
and some examples can be found in [TBJ [EJ E3 HZj- 

The class of Hamiltonians (|4.9|) is invariant with respect to canoniacal transformations 
of the form 

Pi = hpi + hp2 + h, P2 = hpi + hp2 + h, yi = 4>, qi = 4>, 

where hi, ki, <p, <f> are some functions of q±, q 2 . It is easily seen that the functions k±, k 2 , k\, k 2 
can be uniquely expressed through (J), (j) from the condition {p a , qp] = 8 a p\ 

h = ^w-\ k 2 = -^W~\ h = -^w~\ k ™ 



dq 2 dq 1 



where 

W 



dq 2 dqi dqx dq 2 

The fact that the coefficient of p\p 2 in the transformed Hamiltonian ()4.9jl vanishes is 
equivalent to 



Analogously, the condition 



A ^ *%W> + = m *» 1 ^ + ^ 1 (412) 



/2 ^yi uyi t/y2 
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guarantees that the coefficient of ftift 2 in the transformed integral (J4.1(Jj) is equal to zero. 
The existence of a solution (ft, cf> for system (j4.11|) . ()4.12|) is obvious. Our goal is to find the 
functions (/), (f> for the pair ()4.4|) . ()4.5|) explicitly. 

Let us introduce the Kowalewski variables Si(qi,q2) and s 2 (qi, q 2 ) as the roots of the 
equation 

W( qi , q 2 ) s 2 + 2Z(q 1; q 2 ) s + W(q u q 2 ) = 0, 
where W and Z are given by ()4.7|) . (j4.6|) and the polynomial W is defined by f)4.8j) . 

Proposition 3. In the variables s\, s 2 the functions H and K have the form (|1.3JI - (|1.5J) . 

where polynomials S and R are defined by the formulas 

S(x) = — (x + A — ai —a 2 — a 3 )(x + X + ai + a 2 — a 3 )(x + X + ai — a 2 + a 3 )(x + A — «i + a 2 + a 3 ) 
and 

2^ r • • • 

u = -{—-v), S = ji-j 2 , v = 3\+3i- 
4 5 

The values h and k of the integrals H and K are related to the constants e\,e 2 from (jl.2j) 
by 

h = ei + -^aS"(0), k = e 2 + ^aS"\0). 

Thus the general scheme from Section 3 can be applied to the Schottky-Manakov spinning 
top. 

5 The Kowalewski gyrostat. 

In this Section we use some formulas from [21] that describe a dynamical system for the 
Kowalewski gyrostat in the Kowalewski variables. 

The Hamiltonain structure for the gyrostat is defined by the e(3)-Poisson brackets 
{Mi ,Mj} = e ijk M k , {Mi , 7^ } = e ijk ^ k , {7^ , 7j } = 0, 

where Sij k is the totally skew-symmetric tensor. The brackets possess the Casimir functions 

3 3 
A = Y,ll B = Y,lkM k . (5.1) 

k=l k=l 

The Hamiltonian for the gyrostat is as follows: 

H = \^Ml + M\ + 2M 3 2 - 2AM 3 ) + c 7 i , (5.2) 
11 



where c and A are constants. An additional integral of motion is given by 

K = £i6 + 4a((M 3 - \)z lZ2 - + z 2 )c 73 ), (5.3) 



where 



and 



We put 



fi = z\ - 2c(7i + 272) , 6 = ^2 - 2c(7i - 272) 
z x = M l + iM 2 , z 2 = Mi - iM 2 . 
22) = z\z\ - 2h(z 2 + z l) ~ 4cfe(2i + 2 2 ) - 4c 2 a + k. 



Here a, 6, fa and & are values of integrals ()5.1|) . ()5.2|) . and (|5.3|) . The Kowalewski variables 
are defined as follows 

_ R(zi, z 2 ) ± ^(^lj^O-R^j^) 
Sl ' 2 - 2(Z! - , 2 ) 2 • 

It turns out 21 that 



2 VVi ^2/ 2 
5 = (2h + Sl + s 2 )A 2 - \V^~2 (^ + ^-)+(s 1 -s 2 )(^- - s lS2 + h 2 . (5.5) 



4 \ip 1 (p 2 J \ (pi (p 2 

Here ipi = S(s { ), 

S(s) = Ag i -Shs 2 + Ah 2 s-ks + Ac 2 a s + 4c 2 b. 
Substituting the expressions ()2.8|) 

si = -7 + s 2 = -t (v + 1) 

in (|5.4|) . ()5. 5|) for the velocities, we get just ()2.5|) . where 

fc(x) = 4(x + /i) 2 + A5 2 (x - 2h) - k, X = i6. 

Thus in this case the scheme for the separation of variables described in Section 3 is ap- 
plicable. The curve $(y,£) = is of genus 3. The differentials ujx-,w 2 ,u)± form a basis of 
holomorphic differentials. 

Notice that the substitution 5 = —iX in (|3.1j) results the change arcianh — > arctan and 
the action remains to be real. 
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6 The Clebsch spinning top. 



The Clebsch spinning top is defined by the Hamiltonian 

H = \{Jl + Jl + 4) + \{\ix\ + X 2 x\ + X 3 x 2 3 ) 
which commutes with respect to the e(3)-Poisson brackets 

{Jit JjJ i^ijk-Jki Xj} 0, {J'ii %j} i^ijk^k 

with the first integral 

K = (XiJi + A 2 J 2 2 + A 3 J 3 2 ) - A!A 2 A 3 (^ + f + 

M M A3 

Let us fix the values of the Casimir functions as follows 

x\ + x\ + x\ = a 2 , J1X1 + J 2 x 2 + J3X3 = I. 

Using the parameterization 

J i = \pi( l -(li) + \p2iX-ql) + - qi, h = ^-pi(l+qf) + ^-p 2 (l+q 2 )-i- qi, J3 = Piqi+P2q2~-, 
z z (x z z a a 

1 - gi<?2 . 1 + <?i<?2 qi + q2 

Xi = a , x 2 = ta , £3 = a , 

qi - qi " qi- q2 qi - q2 

we express H and K in terms of canonically conjugated variables pi,qi, p 2 ,q 2 - 

1, N9 I . . a 2 W(x,y) 1.. 9 I 2 , 

H = - 2 (x- vfvm + - ft(x - y) + - + -(Aa 2 + -) 



4 v v^i. v*,^. ™ v ,*„ 4a ^, V, 2a r, xv (:T - ?/) 2 12a 2 3d 2 

where 

= (Ai - A 2 )(:r 2 - l) 2 + 4(A 3 - A 2 ):r 2 , 
W(x, y) = A x (x 2 - l)(y 2 - 1) - A 2 (x 2 + l)(y 2 + 1) + 4A 3 xy, 
W(x, y) = A 2 A 3 (x 2 - l){y 2 - 1) - A 3 A!(x 2 + l){y 2 + 1) + AX^xy + «(a; - y) 2 , 
A = Ai + A 2 + A3, k = AiA 2 + A 2 A3 + A3A1. 
We define Kowalewski variables as the roots si(gi, q 2 ) and s 2 (gi, q 2 ) of the equation 

(gi - g 2 ) 2 s 2 + W(q u q 2 ) s + W(gi, q 2 ) = 0. 
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Associating the t-dynamics with K, we get equations (|2.3j) . (|2.4j) for si, s 2 , where 

s-(o = 4(e-A 1 )(e-A 2 )(e-A 3 ). 

Our general procedure for separation of variables from Section 3 leads to curve (jl.ljl . where 

s 6 = 0, 1(0 = ^, 4 = a 2 ^ 2 + (2 ei - a 2 A)£ - e 2 . 

Note that for the Clebsch case we have 6 = i4-. One has to replace arctanh by arctan in 
the formula (|3.1j) and after that the action remains to be real. 



7 Case S = -1 

In the case = — 1 common solutions to the pair (|1.3p . f)l .5|) satisfy the Gibbons- Tsarev 
equation (J2.7|) . In this case the curve has the form 

-Aa5 2 Y 4 + (4ei + 4£e 2 + 4£ 2 a)y 2 + 1 = 0. 

This curve of genus 1 can be parameterized by the Weierstrass function as follows: 

Y = -JL. y {= _^_jl4C_i_ + 1 1 



16a 2 x — (3 3 ' 2a 16 a 2 \x — f3\ x — [3 2 % — @ 3 J ' 

where x = p(z,g 2 ,g 3 ), V = p'(z,g 2 ,g 3 ), 

8ct 2 1 16ct^ 

y 2 = 4(x - (3i)(x - f3 2 )(x - fo), (3 3 = — (4aei - e 2 ,), /?i, 2 = --/% ± — ^— v 73 "- 

The meromorphic integrals in 

3 3 

dt = ^ w l (£n) > rfr = 5^ <4»(£n) 
n=l n=l 

can be taken. The result is given by 

3 3 

* + * + * = «"»■ ^^E'^^f). '--f'-^Efc-Ajk-A) 

where a* = p(z i} g 2 ,g 3 ), yi = p'(z h g 2 , g 3 ). 
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